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Abstract. In the former part of this paper, we give functional equations for Barnes 
multiple zeta-functions and consider some relevant results. In particular, we show that 
Ramanujan's classical formula for the Riemann zeta values can be derived from func- 
tional equations for Barnes zeta-functions. In the latter half part, we generalize some 
evaluation formulas of certain series involving hyperbolic functions in terms of Bernoulli 
polynomials. The original formulas were classically given by Cauchy, Mellin, Ramanujan, 
and later recovered and formulated by Berndt. From our consideration, we give multiple 
versions of these known formulas. 
MSC Numbers: 11M41, 11B68. 



1. Introduction 

Let N be the set of natural numbers, Z the ring of rational integers, Q the field of rational 
numbers, M the field of real numbers, C the field of complex numbers, and No := NU {0}. 
We begin with the classical work of Cauchy |10] who studied the series defined by 
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C 1 - 1 ) y • w \ (seZ), 

z — ' sinh [mir )m s 

m£Z\{0} v ' 

where sinh a; = (e x — e~ x )/2. He showed that several values at s = 4k + 3 (k G No) can be 
written in terms of ir. After his work, this series was considered by Mellin, Ramanujan, 
and several other authors (see [HI El dSl EZ]), and the following fascinating formula was 
proved: 

/, 2) V =(2i) 4W Vf ^+1 -^(1/2) j?4fc+4-2j(l/2) 

1 ' 1 ^ sinh(m7r)m 4fc + 3 1 1 ^ > f 2 j)\ (4fc + 4-2j)! 

m<EZ\{0} V 1 j=0 \ J J J! 

for k G No, where Bj(y) is the jth Bernoulli polynomial defined by 

(see [H]). As a result related to (11. 2p . it is also known that 

/-, 4 x coth(mTr) _ ( 4k+3 ^ , i+1 B 2j (0) B 4k+4 _ 2j (0) 

[ ' L> m«+3 ^ l > ( 2j )\ Uk + 4-2j)\ 

m£Z\{0} j=o JJ 

for k G No, which is written in Ramanujan's notebooks (see Berndt [51 (25.3) p. 293]), 
where cothx = (e x +e~ x )/(e x — e~ x ). In fact, (11. 4p can be easily derived from Ramanujan's 
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famous formula (see Berndt [8, p. 275]): 



a 



1.5) 




l)k 2N 



+i 



_ 9- A \ f -| -^2Af+2-2fc(0) N+l-k ok 

~ 7 ' (2fc)! (2iV + 2-2A;)! P ' 

where iV is any non-zero integer, a and /3 are positive numbers such that a(3 = tc 2 and 
C(s) is the Riemann zeta- function. 

In the 1970's, Berndt [U E] studied generalized Eisenstein series and proved transfor- 
mation formulas for them. Using this result, he gave a family of evaluation formulas for 
certain Dirichlet series in [6j[7], including ( 11. 2p . ( II. 4p and ( II. 5p (see also Remark [6 . 6 [) . 

What is the meaning of the above infinite series involving hyperbolic functions? We can 
find that they are connected with Barnes multiple zeta-functions. In fact, in the former 
half part of this paper, we study functional equations for the Barnes zeta-functions and 
we show two expressions of the Barnes zeta-functions or their residues at integers. We 
observe that Ramanujan's formula (11.51) (and hence (11.41) ) can be deduced by combining 
these two expressions in the double case. Hence in the multiple cases, the combination 
of these expressions may be regarded as generalizations of Ramanujan's formula (see 
Corollary E3D. 

Motivated by this observation, in the latter half part, we first give a very general form of 
evaluation formulas (see Theorem 15. ip . which is out of the frame of Barnes zeta-functions. 
From this form, we deduce a certain explicit evaluation formula with a parameter y G [0, 1] 
(see Theorem 16. 1 p which may be regarded as a relation of several Barnes zeta-functions 
at non-positive integers. This formula especially implies ( 11 .2[) and (ll.4p (see Corollaries 
16.21 and l6.3p and also implies a lot of presumably new formulas, for example, 
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mGZ\{0} 



sinh (rmri / p) 2 m 4 2835 
coth(m7ri/p) 2 62 



7T 4 , 



m 4 2835 



vr 4 , 



where i = \/— T and p = (—1 + y/—3)/2, the cube root of unity, and the same type of 
formulas including higher power roots of unity (see Corollary 16. 4[ Example 16. 5p . 

2. Functional equations for Barnes zeta-functions 

For 9 G R let H(9) = {z = re i{e+ ^ G C | r > 0, -tt/2 < <p < vr/2} be the open 
half plane whose normal vector is e l9 . We recall the Barnes zeta-function defined by the 
following multiple Dirichlet series: 

OO OO j 

(2.1) C»(s, a; wi, • • • , £«>«) = V ■■■ V — ■ ■ -, 

^ ^ (a + wimi H h u n m n ) s 

mi=U m n =0 

where all a, oj\, . . . , u n G H(6) for some 9. Then it is known that this Dirichlet series 
converges absolutely uniformly on any compact subset in 9fts > n. 
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Assume at first that Jis > n. For x G H(9), we have the formula for the gamma 
function 

(2.2) x~ s = —— / e- xt f- l dt. 

r (s) Jo 

Since 

(2.3) a + WiMi + • ■ ■ + u) n m n G H{6) 

for mi, ... , m n G N, we can apply (12. 2p to each term in (12 .ip to get 



oo oo „ e 

Cn(s, a; w x , . . . , u n ) = ^ " " fp\ / 

mi=0 m„=Q ^ ' ^° 



e 



-(a+QJimH hUnmn)^S-l^ 



r (s) 7 (e Wl * - 1) ■ ■ ■ (e"»' - 1) 

1 /" hu; n -a)t 

= - / - t'^dt 

T(s)(e 2 ™ - 1) J c(0) (e"i* - 1) ■ ■ ■ (e<"»' - i) 

where the argument of t is taken in —9 < argt < —9 + 27T and C{9) is a contour which 
starts at e~ ie oo, goes counterclockwise around the origin with sufficiently small radius, 
and ends at e~ l9 oo. Let < yi, . . . , y n < 1 and put 

a = 0(2/1, ...,?/„)= wi(l — 2/1) H h w„(l - y n ) G F(0). 



Then 



C„(s,a(|/i, . . .,y n );u u ...,u n ) 

p (uiyi-\ \-u n y n )t 



(2.5) ?(s){e 2 ™-l) J c{e) (e^-1) 

T-rn —1 

n i= i^ 



T(s)(e 2 ™ - 1) 



n It. 

Jc(ff) j = i 

If t G C(#) is sufficiently far from the origin, then ^R.(ujjt) > (1 < j < n). Therefore the 
integral on the rightmost side converges absolutely uniformly on the whole space C, so 
(12. 5 p gives the meromorphic continuation of ( n (s, d(y); co\, . . . , u n ) to the whole space C. 

In the following, we assume that n > 2 and Q(uj/uk) 7^ for any pair (j,k) with 
j ^ k. From the above integral expression we obtain the following functional equations 
for Barnes zeta-functions. When y± = • ■ ■ = y n = y, we write a(yi, . . . ,y n ) = a(y) for 
brevity. 

Theorem 2.1 (functional equations). We have 

(2.6) C n (s,a(y);wi,...,u) n ) 

2<fti ^-^ /-JL e ( 2m7T iuj/u k )y 



E E ">'(n ^Mn )p™. r 1 )' "« 



r( s )(e 2 ™ - 1) 

V 7V 7 fe=l meZ\{0} 3=1 



where the right-hand side converges absolutely uniformly on the whole space C z/0 < y < 1, 
and on i/ie region 9fts < ify — O. 
In particular, if y = 1/2, we have 

(2.7) Cn(s,(wH hw n )/2;wi,...,w„) 
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12™ n 1 

•2^r( S )(e^-l)£ E ^ 1 (- 1 ) m (n sinh(m7r ^. M) )( 2 ™^ 1 )^ 1 - 

WV 7 fc=l meZUOl .7=1 V ' ?/ 7 



m€Z\{0} j = 



In the case y = of Theorem 12.11 the series expression (I2.6P is valid only for Jis < 0. 
In order to remove this restriction, we decompose the series into the terms involving 
the Riemann zeta-function and the remaining parts. For k G {1, . . . ,n}, let I k = {j G 
{1, . . . , n} \ {k}\Q(uj/co k ) > 0} and I k = {j G {1, . . . , n} \ {k}\$s(uj/u)k) < 0}. Let 



(2.8) 8{J) 

for J C {1, . . . , n}. 
Corollary 2.2. We have 

(2.9) ( n (s,a(0);uJi, . . . ,u n ) 
2iri 



(J^0) 
;-l)«+i (j = 0) 



= " rf^fe^-i^ ^i EC 2 ™**"* 1 )' iQl e2m ^ jM . _ - 

v A ; fe=l Lm>o v j=l 

m>0 ^ j'=l ' 

+ 5(/ fc -)(27rzo; fe - 1 ) s - 1 C(l - a ) + 5(/, + )(-27rio; fc - 1 ) s - 1 C(l - s) 

where the series in the right-hand side converge absolutely uniformly on the whole space 
C. 

Proofs of Theorem 12. II and Corollary 12.21 will be given in Section 4. 
In the following, the empty sum should be understood as 0. 

Corollary 2.3. For I G Z (or / > n ify = 0), we have 

(2.10) J2 E (II ry )^ 1 )^ 1 ^ 

fc=l mgZ\{0} j=l 

_^y-n+i 

Cn(w-^a(y);^i,---,^n) (i>n) 

-(n-Z- l)!Res a=n _jCn(«»a(3/);wi,...,tt; n ) (Z <n). 

On £/ie other hand, this is equal to 

(2.11) - y fr^Mo;^" 1 . 

rai,...,m n >0 j=l J 
miH hm„=l 

Proof. For G Z, the expansion 

rwia^i) = ftg + (.-t) (*<0) 



2vri 



(A;-l)!(s-A;) + 0((s-A;) 2 ) (k > 0) 
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holds. Using this and Theorem 12. 1[ we obtain (I2.10p . On the other hand, by use of the 
integral representation (I2.5p . we see that the left-hand side of (I2.10p is equal to 



(2.13) 



n 

i=i 



Res 

i=0 



J[F(ujt,y)y 



-i-i 



which yields fl2TTPp . 

Similarly we have the following. 
Corollary 2.4. For I £ Z \ {n}, we have 



□ 



k=l 



\m>0 



m>0 



-2m7rzo;i 



-l\n-/-l 



n 

5=1 

n 



,2miviujj /u)). J 



- 5(4") 



(2.14) 



+ 5(/ fc -)(2^ 1 )^- 1 C(l-n + /) 



+5(4 + )(-27r^ fc i r i - i c(i-n+/) 



(Z-n)! 
-fn-Z- 



( n (n - /,a(0); wi, . ..,a; n ) 
l)!Res s= „_iCn(s,a(0);wi, 



(/ > n) 
(I < n) 



Err ^m,(0) m,-i 
11 ^ w 



mi,...,m„>0 j=l 
TOlH hm n =i 



In the next section we will show that Ramanujan's formula ( ll.5p is a consequence of 
the case n = 2 of f)2.14p . Therefore Corollary 12.41 can be regarded as generalizations of 
Ramanujan's formula. 

Here we give historical remarks. A kind of functional equations for the Barnes zeta- 
functions was first proved by Hardy and Littlewood [13] in the case n — 2, and a general- 
ization to the case of general n was discussed in Egami's lecture note [12]. Our proof of 
Theorem 12.11 is essentially the same as those of them. On the other hand, by calculating 
explicitly the residue f!2.13p . we showed an expression of ( n (n — I, a(y); cu\, . . . , co n ) or its 
residues in terms of Bernoulli polynomials. This type of results is also classical, already 
studied by Barnes himself [2], [3] (see also [T], [15] )• In this sense, both of the two equali- 
ties in ( 12.1 Op and (12. lip are classical. The novel point in the present paper is to combine 
these two equalities. A consequence of such a combination is the observation concerning 
Ramanujan's formula in the next section. 



3. Ramanujan's formula 



In this section, we show that Ramanujan's formula (ll.5p can be obtained by combining 
two equalities given in Corollary 12.41 In Corollary 12.41 consider the case n = 2, ui = a 1 / 2 , 
U2 = ifrl 2 with a ,j8eR, Let N e Z \ {0}. 
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The last member of (12.141) is 
(3.1) Q = -(^V^T ^^yf aN)/ 2 ^, 

In particular, for I = 2N + 2, we have 

Cm+2 = -^^Y 1 X>" (2^+2 - j)! B f ] ^ N+2 ^ /2 P /2 
(3 ' 2) nZ 

N+l-k ok 

(2N + 2-2k)\ (2k)l"' 



(j 1/2 01/2^-1 B 2N +2-2k(0) B 2k (0) N+l-k oh 



k=0 



where we have used _B 2 j+i(0) = (for j > 1). 

Next we compute the first member of (I2.14p . Since 5(I±) = <5(i^~) = and <5(/ 1 ~ 
(5(7^) = — 1, this is equal to 



J2(2mma-^r l ( t„ 1/2 _ + l) 



m>0 

+ a -V2^ ( _ 2m7r ^-i/2 ) i-/ 



m>0 



(3.3) -ip-v^^y- 1 

EW")--( r J rw _, + 1) 



m>0 

Q- 1 / 2 (2«a- 1 / 2 ) 1 -'C(/ - 1) + ir 1/2 (-27rr 1/2 ) 1_ 'CC - 1). 



a -v^ (2mxia -i/ 2)1 -._i 

m>0 



In particular, in the case a 1 / 2 /? 1 / 2 = ir, we see that this is equal to 

, 1 

e 2m/3 _ ]_ 

+ a -v^ ( _ 2m , iQ -V 2) i-«_i_ 

m>0 

_ ir v 2 £ (2m . r v T -._L_ 

m>0 

(3.4) +ir i/^(_ 2m , r i/ 2 )-_L_ 

m>0 

- a" 1 / 2 (2 7 rza- 1 / 2 ) 1 -'C(/ - 1) + ip- 1/2 (-W 1,2 ) 1 - l C(l ~ 1) 

oo 

= (2.)-¥-' a <'- 2 )/ 2 (-C(/ - 1) + ((-l)'- 1 - 1) £ _ !) 



1 



m' 1 



- (^^^((-IJ'CCI - 1) + (1 - (-D 1 -) £ (e^-lW- 

m=l ^ ' 

Further in the case I = 2N + 2, we see that (13. 4p reduces to 
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(3.5) (m)-\2n)-™(-af(-\a2N + 1) - £ 



3 2m/3 _ l) m 2iV+l 



m=l 



oo 

+ ( 7 rz)- 1 (2 7 r)- 2 ^(^(2iV+l) + ^ 



3 2mo 



— l)m 



27V+1 



)- 



By equating (13. 2p and (13.51) . we finally obtain 
(3.6) -2 , 9Ar . 9 _ 9M , , 9n , 



iV+l-fcofc 



fc=0 



(2iV + 2-2A;)! (2fc)! 
-/3)-^(--C(2iV + l)-^ 



0* 



3 2m/3 



m=l 



l)m 



2AM-1 



oo 



' (e 2ma - l)m 2JV + 1 



which recovers (jl .5j) . 



4. Proofs 



Proof of Theorem \2.1[ For z G C and e > 0, let D(z,e) be the closed disk whose center 
is 2 with radius e. We use the notation x + = max{i, 0} for x £ R. We first note that 
for any e > 0, there exists M' = M'(e) > such that for t G C \ Umez D(2miri, e) the 
inequality 



(4.1) 



e* - 1 



holds. Hence for t G C \ Uf=i Umez D{2nmibJ- , e), we have 



(4.2) 



j[]F(^f.//) < M|t| n eS"=i^*)»-»(« J -0 + ) 



with a certain M = M(e, uii, . . . , u n ) > 0. Since 

(4.3) R{ujt)y - U(ujt) + < 

for 1 < j < n, we see that there exists T = T(e) > such that for all t G C with \t\ = 1, 

(4.4) - (^)+) < " T - 

3=1 

Hence we see that for alH G C \ Uj=i Umez D(2mniuJ 1 , e), 

n 

(4.5) |jjF(w^,y) < M|t| n e- T| *L 

3=1 

If < y < 1, then we can choose T > 0. In fact, since ^(uj/uk) 7^ for j ^ k and n > 2, 
for any £ with \t\ = 1 we find at least one j for which di(ujjt) 7^ holds. Then, using 
< y < 1 we see that $t(tOjt)y — $l(ujt)+ < 0, from which T > easily follows. 
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From (14. 5p . we see that the integral on the rightmost side of (12.51) converges to when 
the radius of the contour goes to infinity if < y < 1 or, y = with 3fts < 0. Namely, 
there is a sequence Ri — > oo such that 



t s - n - 1 ||dt| = 0. 



n 11 

(4.6) lim / iTTF^y) 

Hence we can calculate the integral by counting all the residues on the whole space. 
Since by the assumption the poles of the integrand is all simple except the origin, we 
obtain 

(4.7) C n (s,a(y);wi,...,u) n ) 

1 r P (^iH \-u n )yt 

= - / - f^dt 

T(s)(e^-l)J c{e) ( e ^-l)...( e ^-l) 



2 " T) ± e .in^wv-, 



r( 8 )(e Ma - 1) ^ 

v yv ' k=l m£Z\{0} j=i 

whose absolute and uniform convergence follows from the explicit form of the series. 
Therefore we obtain (12. 6p . □ 

Proof of Corollary \2.B. We observe that as m — >■ +oo, F(2irimujj/ujk, 0) = 0{m) if j G J^" 
while F(2iiimuj/cok,0) decays exponentially if j G I k . Thus we see that if I k ^ 0, the 
series 



1 

£(n 



nmiuj k 



m>0 j=l 

converges absolutely uniformly for whole s G C. 
Next consider the case I k = 0. We have 



( 4 - 8 ) Edl ^-M,! )^"^ 1 ) 



m>0 j=l 



E(- i ) w " i ( 2 ^" i ) s ~ i (ii( i + e -w.-M _ i 

m>0 i=l 

^(-ir- 1 (2m 7 r^ fc 1 )- 1 [l+ £ (J] 
(-ir 1 (2 7 rza; fe - 1 )- 1 C(l- S ) 

+(-ir 1 E E^^^^dl 



m>0 \ Jc{l,...,n}\{fc} j£J 

|J|>1 



1 



g—lmiriujj jui^ ^ J ' 

,/C{l,...,n}\{fc} m>0 jeJ 
|J|>1 



Since all j G J^", we see that the rightmost side of the above can be continued to the 
whole of C, and is equal to 

(4.9) (-ir- 1 (2vr^T 1 ) s - 1 C(l- S ) 
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+ i> m ^'r'((n e2m j llm _ j - (-1)-). 

m>0 ^ j=l ' 

For the series with m < 0, by exchanging the roles of l£ and I^T, we have the same 
type of conclusions as follows: If l£ ^ 0, the series corresponding to m < converges 
absolutely uniformly for whole s G C, and if l£ — 0, 

n -. 

( 4 - 10 ) E (II ^v^i ) (a^,- 1 )'- 1 

m<0 3=1 

m>0 V j=l 7 

□ 



5. A GENERAL FORMULATION 

In the previous sections, we established some relations between Barnes zeta-functions 
and certain series involving hyperbolic functions (see (12 .7p ). In order to study this re- 
lationship further, it is convenient to introduce a general framework to evaluate more 
general series. 

Let g(t) be a meromorphic function on C which has possible poles only on 1-niL. For 
example, we will consider g(t) = it/2)/ sinh(t/2) (see ( I6.8P ). 

Let ?i6N with n > 2 and r) = e m ^ n , that is, the primitive 2n-th root of unity. Let 



G(*)=n^)- 

We assume that there exist real numbers 71, 72 with 71 > and a small positive number 
e such that 



\G(t)\ < 7 i|t| 



for all t G C \ ljr=o 1 UmGZ D^rrntirf , e). Then we have the following theorem, which is a 
simple consequence of residue calculus, but is a key result in the present paper. 

Theorem 5.1. For h 6 Z mft /i + 72 > 1, 

/ x \ 

n— 1 n 



(5.1) 



m£Z\{0} 



V 



«=0 



Res .#(-£) 

£=2m7Tt 



/ 
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In particular when g is an even function, 



fn-\ 



Z h V (2mni)- h [T\g(2mmif) ) Res ,g(t) 
^ ^ mez\{o} \j=i / 

where 

n _i In if h = 1 (mod 2n), 

(5.3) Z h = ^r)^ 1 -^ = < ifh^l (mod 2n) and 2 J(h, 

j=o [ j^tk ifh^l (mod In) and 2 | h. 

Proof Let h G Z with h + 72 > 1. For i? e N, we have 
(5.4) 

ResjcmrH = — [ G(t)r h dt 



= V V Res + — / G(t)r h dt 

t^o<H<^ =2mW W| t | =27r * + 2 £ 
Since G(t)t _fc < 7 1 |t|- /l ~ 72 , we have 

/ G{t)r h dt < [ \G{t)r h \\dt\ 

J\t\=2wR+2e J\t\=2irR+2e 



(5.5) <7i/ 1*1-^1*1 

J|t|=27ri?+2 £ 

< 27T7 1 (27r J R + 2e)^ 72+1 
->■ 

as R — > 00. Hence by letting R — )■ 00, we obtain 

n— 1 

(5.6) Res{G(t)r h } = -J2v l(1 - h) £ WO* - *} 

mGZ\{0} 

because 



(5-7) < 71 / Itl-^ldtl 



=0 






1 


< 






27 


< 


7i 




2tt 



Res G( V l t)t~ h <±-f \G{rft)r h \\dt\ 



\t—2miri\=2e 



=2e 



< 2£7 1 (27r|m| - 2e)^ 72 
and hence the convergence is absolute. 
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Since rj is the primitive 2n-th root of unity, we see that for I 6 Z with < I < n — 1, 
the residue of 

n 

t- h G( V l t) = t~ h Y[ g{r] j+l t) 
j=i 

at t = 277111% is equal to 

n 

(2rmri)- h TT g(27mri7] j+l ) x Res ,g(-t), 

3=1 
j+l^n 

which gives (15. ip . 

In particular when g is even, we have G(t]H) = G(t) (0 < I < n - 1), because #(r/ r t) = 
g{—Tf~ n t) = g{i] r ~ n t) f° r ^ + 1 < r < 2n. Therefore 

n-1 

E^ E t toa.{T h Gtft)} = Z h E Res.{r h L7(t)}. 

' * ^— » t=2rnirt ' * t=2rmri 

1=0 mGZ\{0} mGZ\{0} 

This completes the proof. □ 

The same idea as in the above proof is also used in [Hj in a different situation. It is to 
be noted that the original method of Cauchy [10] is essentially similar. 

6. Explicit formulas 
We recall the Bernoulli polynomials {Bj(y)} defined by (11. 3ft : 

For t/GR, let 

„, x F(t,y) + F(-t,y) t e%-V2) + e -%-i/2) 

rfin = 2 = 2 - e-/ 2 

1 j t cosh(% - 1/2)) 

~ 2 sinh(t/2) ' 

Then we see that 

(6.2) H(t,y) = Y< B ^y) 



(2m)!' 

m=0 v y 

It follows from (16. ip that H(t,y) has simple poles at t = 2m7ri (m 6 Z \ {0}) and its 
residue is 

2m7ri e 2 ""™^- 1 / 2 ) _|_ e - 2 ««r%-i/2) 
(6.3) Res H(t,y) = -, r = 2rrmi cos(2mny) . 

V ' t=2mni V 2 (-l) m V ^ 

By (15. 2p . we have the following result which includes the known formulas (11.21) and (11.41) 
given by Cauchy, Mellin, Ramanujan and so on (see Section [1]). 

Theorem 6.1. Assume < y < 1 and pGZ, or y = 0,1 and p > n/2. For nGN with 
n > 2 and r\ = e n ^ n , 

(6 A) Z V cos(2m7iy) ( cosh (2m?r irf (y - 1/2)) 
10- ) 2 P +i 2^ m 2 P +i~n I 11 sinh(mW) 

mGZ\{0} \j=l V ' ; 
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_ 2"- 1 (27T»)fr+ 1 -" v " B 2mv (y) {v _ l)mu 

n(n-l)/2 11 ( 2m;/ )! ' 

mi,...,m n >0 u=l v ' 
toiH \-m n =p 

where is defined by (15. 3p . Furthermore, assume < y < 1 and p > n/2. Then the 
both sides of (16. 4p «s a/so egna/ to 

/o •\2p+l-7i/_i \l-n n_1 

(M ( l 2 „_ ( w) , E E Cn(n-2 P , $> % ,1,, ,...,,«-). 

2/oe{y,i-2/} Vn-ie{v,i-v} i=° 

Proof. Since H(—t,y) = H(t,y), we can apply (15. 2p with g(t) = H(t,y) and h = 2p + 1. 
In this case 72 is arbitrarily large for < y < 1 and 72 = — n for y = 0,1, and hence the 
condition h + 72 > 1 is satisfied because we assume p>n/2ify = 0,l. By ( 16. 2p . we have 



(2m)! 

m=0 v 7 



Hence we obtain 

n-l 



Res ^ r*- 1 J] Htft, y)\= E f[ ^SV" 



(f-l)m» 



i=0 J mi,...,m B >0 
mi H hm n =p 



■~ (2m,)! 

_+---+77l„=J 

Therefore, by using (16. ip and ( I6.3p . we have 



m6Z\{0} \j=l 



-f -B2m„(j/) ^2^-l)?nv 



V TT 

^ 11 (2m„)! ' 

mi,...,m n >0 ^=1 



il,...,m n: 

7711 H h77ln=P 



Thus we obtain (EU). 

Assume < y < 1 and p > n/2. Then /i = 2p+l>n+l. Note that 

e 1 — 1 e* — 1 



Then we have 



(6.6) 



i=o j=o 

n-l 

= 2 ~ n e ••• e n /-(//-/-//,)- 

yoG{y,l-j/} j/ n _iG{y,l-j/} i=0 



Hence 

Ra»{G(t)r fc } 
(6-7) 2 _ 

yo&{yA-y} y n -i£{y,i-v} |t|=e i=° 



O — 71 x» 71 1 

s E - E / (n^w))***"-^ 1 * 
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n-1 



Since < y < 1, we see that 1, rj, . . . , rf 1-1 and X]j=o ^'(l ~~ Vj) are belonging to the half 
plane H(9 n ), where 9 Tl = n/2 — 7r/(2n). Therefore, deforming the path |i| = e to C(8 n ), we 
find that each integral on the right-hand side of (16. 7p is of the same form as the integral 
on the right-hand side of (12. 5p with s = n + 1 — h. Using (12.51) and (I2.12p . we obtain that 
the right-hand side of (16.71) is equal to 



2~ n„n(n— l)/2 

2tH 



E - E 

yo£{y,i-y} y n -ie{y,l-y} 



QW — 1 \h-n-l ™-l 

(h-n-l)\ Cn{n + 1 - ^ E ^ - !' * " " " > ^ 

<2^— n / j \ n— l^n(n— 1)/2 



(h-n- 1) 



n-1 

n-l\ 



W)G{y,i-y} j/n-ie{y,i-»} i=o 
which implies ( 16. 5p . □ 

In particular when y = ^ in (16. 4p and (16. 5p . we have the following formula, which can 
be regarded clS cL multiple generalization of ( II. 2p . 

Corollary 6.2. For p E Z, n e N n > 2 and i] = e 7 ™/™, 

(_l)m 



^2p+l E 



(6.8) 



2p+l— n 



2 n-l (27r , )2p+ l-n » ^(l /2) ^ 

n(n-l)/2 11 (2m v )\ ' 

mi,...,m„>0 v=l v ' 



-il,...,m nt 
miH hm„=p 



Furthermore if p > t/ien ( 16. 8 p is egwa/ to 

(6.9) 1 J 2 (2p -n)! C> " ^ 1/(1 " ^ V > " " " ' ^ } " 

Now we consider the case n = 2, r] = e 7 ™/ 2 = z and p = 2k + 2 in (16. 4p . Then we obtain 
the following. 

Corollary 6.3. For k e~N and < y < 1, 

\ cos(2rmry) /cosh (2rmr(y — 1/2)) \ 

m£Z\{0} v \ ) / 

( 6 - 10 ) 



_ / _\4Jfc+3 V s / lV? -+1 ^2j(y) B 4k+4 _2j(y) 

~ { ] } m (4Jb + 4-2j)!- 

In particular when y = ane? y = | ; we obtain (ll.4p and (II. 2p . respectively. 
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Next we consider the case n = 3. Let n = e m ^ 3 = —p 2 with p = e 2m ^ 3 and p = 3k + 3. 
Then we have r/ 3 ( 3-1 )/ 2 = — 1 and -Zgffc+il+i = 3 by (15. 3p . From (16 .4p . we have 

cos(2m7n/) /cosh (2mrri(—p 2 )(y — 1/2)) cosh (2mnip(y — 1/2)) \ 
m 6fc+4 \ sinh (m7rz(— p 2 )) sinh (unrip) J 



(6 - U) 2 2 ( 27 rz)^+ 4 v 

-3 ^ J-- 1 - (2m,)! P 

mi,m2,m3>0 j/=1 v ' 
mi+m2+m3=3fc+3 

Note that p 2 = l/p and p = — l/p — 1. Hence, by using 

sinh(m7rip) = — (— l) m sinh(m7ri/p), 
we can rewrite (16. lip as follows. 

Corollary 6.4. For fceN and < y < 1, 

(_l) m cos(2mny) /cosh ((2miri/ p)(y — 1/2)) cosh (2muip(y — 1/2)) \ 
-6fc+4 ^ sinh (rnvri/p) 2 / 



m u 

m£Z\{0} 



E (2m!)! (2m 2 )! (2m 3 )! P 



(6-12) 4( _ 1)fe(27r)6fc+ 4 B2mM B2mM B2mM 

3 ^ 

mi,m2,r«3>0 
mi+m2+m3=3fc+3 

In particular when y = anc? y = ^, £/ie following equations hold: 
coth(m7rz/p) 2 



(6.13) £ 

m£Z\{0} 



m 6fc+4 



4(-l) fc (27r) 6fc + 4 ^ B 2mi (0) ^ 2m2 (0) ff 2m3 (0) 2W3 
3 ^ >n (2m x )! (2m 2 )! (2m 3 )! 9 

mi ,1712 ,m3>0 
mi+m,2+m3=3fc+3 



(6.i4) ^ 1 



sinh(m7T2/p) 2 m 6fc+4 

m£Z\{0} v ' F ' 

_ 4(-l) fc (27r) 6fc + 4 ^ B 2mi (l/2) i? 2m2 (l/2) iWl/2) m2+2m3 

(2m!)! (2m 2 )! (2m 3 )! P 



mi,m2,m3>0 
mi+m,2+m3=3fc+3 



Example 6.5. From (16.131) and (I6.14p . we obtain 

coth(mvri/p) 2 62 4 

( 6 - 15 ) 2^ ^4 = 



m 4 2835 

m^O 



(6-16) £ 



coth(mvii/p) 2 40247 



m 10 1915538625 

m^O 



7T 



10 



' » sinh mir' 



sinh(m7ri/p) 2 m 4 2835 



1 ^ 



(6.18) - - ' - - - 7 -° 3 ._- Vr 10 . 



sinh(m7rz/p) 2 m 10 1915538625 
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Additionally, by setting (p,n) = (4,4), (5,5) in equation f|6.8|) . we obtain 



(6.19) 



-l r i 5 

7T 5 , 



, Q sinh(m7r) sinh(m7r^8) sinh(m7u£ 8 1 )m 5 37800 



(6.20) 



-l) m 1 



, q sinh(m7r^ 5 ) sinh(ra7nCf) sinh(TrariCf) sinh(m7rz£|)m 6 467775 



7T 



where Cfc = e 2wi/k (k e N). 

Remark 6.6. By using the same method as introduced in this paper, we can recover the 
known formulas, for example, 

(— l) m 



E 

m=0 



cosh((2m + l)7r/2) ((2m + l)/2) 



4fc+l 



2 A- 



E 

m=0 



(27r )4 fc+ i - ff 2j (l/2)E 4fc _ 2j (l/2) 
8 1 (2j)\ (Ak-2j)\ 



cosh((2m + l)V37r/2)((2m + l)^) 6 ^ 1 



_ 1)fc+ i (27r)6fc+ i ™ ^. +i(0) ^^.(Q) /(2j + l)7r 

_ / . , , _ _ _ \ ■ cos 

^ 



,_ (2j + 1)!(6A;-2j)! 

for fc 6 N (see Watson [18] and Berndt [7]), where {i?„(x)} are the Euler polynomials 
defined by 

2e xt ^ „ , t n 

= \ EJx) — . 

e* + l ^ y ! n\ 

More generally, we can give relevant analogues of these results like those in Example 16.51 
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